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(SI ■ ABSTRACT 



Aims. We study the circumstellar environment of the M-type AGB star RT Vir using mid-infrared high spatial resolution observations 
from the ESO-VLTI focal instrument MIDI. The aim of this study is to provide observational constraints on theoretical prediction 
' that the winds of M-type AGB objects can be driven by photon scattering on iron-free silicate grains located in the close environment 

• ' (about 2 to 3 stellar radii) of the star. 

i-J^ I Methods. We interpreted spectro-interferometric data, first using wavelength-dependent geometric models. We then used a self- 

consistent dynamic model atmosphere containing a time-dependent description of grain growth for pure forsterite dust particles to 
reproduce the photometric, spectrometric, and interferometric measurements of RT Vir. Since the hydrodynamic computation needs 
' stellar parameters as input, a considerable effort was first made to determine these parameters. 

-^-^ I Results. MIDI differential phases reveal the presence of an asymmetry in the stellar vicinity. Results from the geometrical modeling 

^ . give us clues to the presence of aluminum and silicate dust in the close circumstellar environment (< 5 stellar radii). Comparison 

between spectro-interferometric data and a self-consistent dust-driven wind model reveals that silicate dust has to be present in the 
region between 2 to 3 stellar radii to reproduce the 59 and 63 m baseline visibility measurements around 9.8 /jm. This gives additional 
observational evidence in favor of winds driven by photon scattering on iron-free silicate grains located in the close vicinity of an M- 
^ , type star. However, other sources of opacity are clearly missing to reproduce the 10-13 yum visibility measurements for all baselines. 

(SI ' Conclusions. This study is a first attempt to understand the wind mechanism of M-type AGB stars by comparing photometric, 

' spectrometric, and interferometric measurements with state-of-the-art, self-consistent dust-driven wind models. The agreement of the 

00 I dynamic model atmosphere with interferometric measurements in the 8-10 fim spectral region gives additional observational evidence 

l/^ , that the winds of M-type stars can be driven by photon scattering on iron-free silicate grains. Finally, a larger statistical study and 

progress in advanced self-consistent 3D modeling are still required to solve the remaining problems. 

c 2 ^ ■ Key words, techniques: interferometric - techniques: high angular resolution - stars: AGB and post-AGB - stars: atmosphere - stars: 

CO ' circumstellar matter - stars: mass-loss 

1. Introduction an outflow. As a response to these findings, Hofner (2008) 

/\, . proposed an alternative solution, involving thermodynamic 

. The mass loss of AGB stars is commonly ascnbed to a two- conditions where iron-free silicate grains are able to reach sizes 

_C^j step process, where the pulsation of the star creates atmospheric ^^e range of O.l-l.O^^m in the close circumstellar environment 

shock waves that push the matter to conditions where both tem- (2 to 3 stellar radii). State-of-the-art radiation-hydrodynamic 

perature and density allow the formation of dust grains. When ^^^^^^ demonstrate that forsterite particles (Mg2Si04) in this 

the dust opacity is large enough, grains receive enough momen- ^j^^ ^^^g^ ^^^^ ^^^^ ^^^^^^ radiative scattering cross-sections 

tum through radiative pressure to be accelerated away from the ^^^^^ ^„ ^^^^^^ (H^f^c^^ 2008) and that the resulting visual 

star and to drag the gas along by coUisions, causing a slow out- ^^^^.j^ colors agree with observations (Bladh et al., 201 1, 

2012b). Recently, aperture-masking polarimetric interferometry 

Although this scenario has been successfuUy tested for C- ^jth jhe NAOS-CONICA VLT focal instrument at near-IR 

type AGB stars by comparing detailed models with observations wavelengths has revealed the presence of scattered flux resulting 

(e.g. Winters et al. 2000; Gautschy-Loidl et al. 2004; Nowotny f^om 0.3//m sized silicate grains in the close environment 

et al. 2010, 201 1 ; Sacuto et al. 201 1), some fundamental open (^^out 2 stellar radii) of three low mass-loss, M-type AGB stars 

questions remain for the case of M-type AGB stars. As shown (Morris et al., 2012), giving strong observational support to the 

by Woitke (2006), silicate grains have to be basically iron-free theoretical prediction by Hofner (2008). 
to form and survive in the close vicinity of the star, resulting 
in insufficient absorption cross sections in the near-IR to drive 

Because of the poor angular resolution of the IRAS, ISO, 

* Based on observations made with the Very Large Telescope and SPITZER satellites (larger than 2"), the mid-infrared excess 

Interferometer at Paranal Observatory under programs 083.D-0234 and of AGB stars observed with these instruments is dominated by 

086.D-0737 (Open Time Observations). the thermal emission from the extended dusty circumstellar re- 
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gions (above 100 stellar radii) where iron-rich silicate dust par- 
ticles (Mg-Fe-silicates in the form of oUvine and pyroxene) are 
thermally stable (e.g. Woitke 2006). A clear way to disentan- 
gle iron-free and iron-rich silicate dust emission is to spatially 
resolve the innermost parts of the stellar envelope. High angu- 
lar resolution observations in the mid-infrared are very attractive 
for this purpose, as they probe the close circumstellar region in- 
side five stellar radii (e.g. Ohnaka 2004), where only iron-free 
silicate particles can form and resist evaporation due to radiative 
heating (Woitke, 2006). By modeling the lOyum interferomet- 
ric emission with self-consistent dynamic atmosphere and wind 
models, we test the theoretical predictions of Hofner (2008) and 
give additional observational constraints regarding winds driven 
by photon scattering on iron-free silicate grains. 

RT Vir is an M-type (Sp.Type M81II), semi-regular variable 
star and one of the brightest water maser sources (e.g. Imai 
et al. 2003; Richards et al. 2011). The mass-loss rate of the 
star, found from water maser emission and CO line emission, 
is around 1 to 5x10"^ Moyr"' (Bowers & Johnston, 1994; 
Olofsson et al., 2002; Imai et al., 2003). Measurement values 
of the gas velocity vary from 6.2kms~^ (Gonzalez Delgado 
et al., 2003) in the SiO emitting region (~ 400 AU away from 
the star) to 7.8 km s"' (Olofsson et al., 2002) in the outer CO 
emitting region. Gas velocities up to 11 kms"' (Richards et al., 
2011) have been measured in the water maser emitting region 
(~5 to 25 AU away from the star). The mid-infrared excess 
observed between 9yum and 13yum in the ISO/SWS data (see 
Fig. 3) indicates that the star is surrounded by a dust shell. 
The pulsation period of the star has been estimated to 155 
days by Kholopov et al. (1985) and reevaluated to 375 days 
by Imai et al. (1997) using the AAVSO visual lightcurve with 
magnitudes varying from V-l.l to 9.7. The Hipparcos distance 
to the star is 135+15pc. The period-luminosity relation of 
Whitelock et al. (2008) gives a distance of 226 pc, in very good 
agreement with the distance of 220 pc derived by Imai et al. 
(2003) using both the statistical parallax, and model-fitting 
methods for the maser kinematics. We assume in the following 
that the distance of 220 pc determined by Imai et al. (2003) 
is the most reUable, since the Hipparcos distance could suffer 
from a systematic bias due to photocenter displacement of the 
star (Chiavassa et al., 2011). The effective temperature of the 
star is estimated to 2900 K (see Sect. 4.2), while the stellar 
diameter has been evaluated to 12.4 mas in the K and K' bands 
(Monnier et al., 2004; Mennesson et al., 2002) and to 16.2 mas 
in the L' band (Chagnon et al., 2002; Mennesson et al., 2002). 
Assuming that the uniform disk (UD) diameter in the K band is 
a good approximation of the continuum photospheric diameter 
(Scholz & Takeda, 1987), and that the distance is 220 pc and 
the effective temperature is 2900 K, the luminosity of RT Vir is 
about 5500 Lq. 



The MIDI mid-infrared spectro-interferometric observations 
of RT Vir are presented in Sect. 2. Section 3 gives a first mor- 
phological interpretation of the object, based on the N-band in- 
terfero metric data. The determination of the most appropriate 
dynamic model atmosphere for RT Vir is described in Sect. 4. In 
Sect. 5, we aim at constraining the dust chemistry in the inner 
wind region of RTVir, in order to test the theoretical predic- 
tions of Hofner (2008). In Sect. 6 an alternative geometry is ex- 
plored. Finally, the conclusions and future perspectives are given 
in Sect. 7. 



2. Spectro-interferometric observations 

The Very Large Telescope Interferometer (VLTI) of ESO's 
Paranal Observatory was used with MIDI, the mid-infrared 
(A-S-13 jum) interferometric recombiner (Leinert et al., 2003). 
MIDI combines the light of two telescopes and provides single- 
dish acquisition images, calibrated spectra, visibilities, and dif- 
ferential phases in the N-band atmospheric window. All the ob- 
servations, with the exception of observation #7 (see Table 1) 
made in HIGH-SENS mode, were performed in the SCI-PHOT 
mode of MIDI, meaning that photometry and interferometry 
were recorded simultaneously. The GRISM mode was used, pro- 
viding a spectral resolution of about 230. 

The observations of RT Vir were conducted in 2009 and 201 1 
with the VLT auxiliary telescopes (ATs) GO-HO, DO-HO, AO-Gl, 
and AO-KO. These configurations provide projected baselines 
ranging from 128 m down to 30 m, probing spatial scales slightly 
above the stellar photosphere (~1.5 stellar radii) to regions lo- 
cated above the dust condensation zone (~6.5 stellar radii). 

2. 1 . MIDI visibilities and differential phases 

The data reduction software package MIA+EWS' was used to 
calibrate the visibility data. We checked that the resulting cal- 
ibrated visibiUties difl'ered less than 5% to 10% between the 
reductions performed with MIA and with EWS. Because EWS 
also provides diff'erential phases, the result of EWS is presented 
in the remainder of the paper. Various instrumental and atmo- 
spheric effects can affect the fringe contrast of the science tar- 
get. CaUbrator stars with known characteristics and dimensions 
must be observed close in time and space from the science tar- 
get to derive the response of the instrument, and properly cal- 
ibrate the interferometric measurements (e.g. Cruzalebes et al. 
2010). During the first period (P83), we obtained measurements 
of the calibrator HD 120323 (F,2^™=255.4 Jy and 0=9.16+0.07 
mas)^ before or after the measurements of RT Vir, while in the 
second observing period (P86), measurements of the same cali- 
brator were obtained before and after RT Vir (see Table 1). The 
data quality has been checked using the MIA Graphical User 
Interface (e.g. Sacuto et al. 2011; Klotz et al. 2012a). Among 
the MIDI measurements, data set #4 shows very broad dispersed 
power spectral densities (PSD) and fringe histogram (FH), prob- 
ably caused by the bad atmospheric conditions during the night 
(seeing>2"). If the PSD and FH are very broad, a very significant 
fraction of the fringe power is distributed outside the integration 
range, leading to a systematic underestimation of visibihty (e.g. 
Sacuto et al. 2011). The corresponding visibihty measurements 
were therefore discarded. 

For P83, the calibrated visibilities from 2009 March 18 (data 
sets #1 and #2) were averaged over both bracketing calibra- 
tor observations obtained during that night, while error bars 
were derived from the standard deviation of the visibilities. 
This was not possible for the 2009 June 03 data (data set #3), 
since no other calibrators were observed in the same mode 
(GRISM+SCI-PHOT) during that night, and a standard mul- 
tiphcative error of 10% was arbitrarily assumed for the visi- 
bilities (Chesneau, 2007). The calibrated visibilities from P86 
were averaged over both bracketing calibrator observations ob- 
tained during each night, while error bars were derived from the 



' http : //www . mpia-hd . mpg . de/MIDISOFT/ , 
http : / /www . strw . leidenuniv . nl/~nevec/MIDI/ 

^ IRAS 12//m flux and angular diameter from the ESO/MIDI cali- 
brator database, respectively 
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Table 1. Journal of the MIDI Auxihary Telescopes observations of RT Vir with its main calibrator. 



# 


Star 


UT date & Time 


0v (Cycle) 


Config. 


Baseline length [m] 


PA [deg] 


Seeing [arcsec] 


1 


RTVir 


2009-03-18 05:03:28 


0.15 (1) 


DO-HO 


59 


74 


1.7 


Cal 


HD 120323 


2009-03-18 05:23:40 




- 






1.5 


2 


RTVir 


2009-03-18 06:03:49 


0.15 (1) 


DO-HO 


63 


73 


1.0 


Cal 


HD 120323 


2009-03-18 06:26:04 




- 




... 


0.8 


3 


RTVir 


2009-06-03 01:33:04 


0.35 (1) 


GO-HO 


32 


73 


1.0 


Cal 


HD 120323 


2009-06-03 02:00:57 




- 






1.1 


4+ 


RTVir 


2009-07-01 00:26:23 


0.43 (1) 


AO-KO 


126 


72 


2.4 


Cal 


HD 120323 


2009-07-01 00:03:48 




- 






2.0 


5 


RTVir 


2011-02-08 08:29:37 


0.00 (3) 


AO-Gl 


89 


113 


0.9 


Call 


HD 120323 


2011-02-08 08:06:59 










0.9 


Cal2 


HD 120323 


2011-02-08 08:55:11 










1.0 


6 


RTVir 


2011-03-06 06:10:26 


0.06 (3) 


GO-HO 


30 


74 


1.2 


Call 


HD 120323 


2011-03-06 05:42:11 










1.2 


Cal2 


HD 120323 


2011-03-06 06:35:04 










1.2 


7 


RTVir 


2011-03-13 07:06:34 


0.08 (3) 


AO-KO 


128 


73 


1.4 


Call 


HD 120323 


2011-03-13 06:46:31 










1.3 


Cal2 


HD 120323 


2011-03-13 07:26:59 










1.2 



Note 1: Calibrators used to calibrate the visibilities are given after the science target. The visual phase 0v of the corresponding observing cycle is 
calculated from Eq. 1. The configuration, the length and position angle of the projected baseline, and the seeing value are also indicated. 
Note 2: ^ means that the data were not exploitable (see text). 
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Fig. 1. Spectrally-dispersed MIDI calibrated visibilities of RT Vir ob- 
served over 5 nights between March 2009 and March 201 1 (see Table 1) 
at baselines 30(green lines), 60(red lines), 90(blue line), and 130 m (pur- 
ple line) from top to bottom, respectively. 



Fig. 2. Spectrally-dispersed MIDI calibrated differential phases of 
RTVir (similar color-code as in Fig. 1) observed over 5 nights between 
March 2009 and March 2011 (see Tab 1). The bold purple line repre- 
sents the largest 128 m baseline measurements (data set #7). 



standard deviation of the visibilities. Figure 1 shows the MIDI- 
calibrated visibilities of RT Vir observed over the five nights (see 
Table 1). 

Error bars for the calibrated differential phases were derived 
in the same way as for the visibilities. Because the humidity and 
seeing levels were stable during 2009 June 03, a typical error 
of ±5° was used for the data set #3 differential phases (Ohnaka 
et al., 2008). Figure 2 shows the MIDI calibrated differential 
phases of RTVir observed over the five nights (see Table 1). 
It is clear from this figure that the 128 m baseline differential 
phases show a significant deviation from zero in the 8-9 jum 
band^, greater than the expected upper-limit uncertainty of 10°. 



^ Differential phase measurements for wavelengths longer than 
12.0/jm are more uncertain owing to the increasing effect of the at- 
mospheric water vapor and edge band noise. Therefore, the differential 
phases larger than the upper-limit uncertainty seen between 12.0 and 
13.0 fim for data set #5 are probably not real. 



A similar differential phase signature was found by Paladini et 
al. (2012) in the close environment of the C-type star R For 

2.2. MIDI and ISO/SWS spectra 

MIDI spectra are derived from the average of the calibrated 
spectra close in visual phase within a given observing cycle, 
while the error bars are determined from the standard deviation 
of the calibrated spectra. Since the calibrated spectrum of data 
set #7 leads to an unusual spectral signature that is different 
from the other calibrated spectra, we decided to discard this 
measurement together with data set #4 (see Sect. 2. 1). 

In the following, the visual phase of the star is estimated 
from its AAVSO lightcurve (Henden, 2009) according to 

(f-rQ)modP 

H = p . (1) 
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^ n n F ■ 
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wavelength (//m) 

Fig. 3. Comparison of the ISO/SWS spectrum of RTVir (black thick 
line) and the MIDI calibrated spectra of RT Vir at phase i^v =0. 1 5 (Cycle 
1) (blue solid line + error bars), phase 0v=O.35 (Cycle 1) (red dashed- 
dotted line), and at phase 0v~O.O (Cycle 3) (green dashed line + eiTor 
bars). 

where t is the observing time in Julian date, Tq-2 454 854 is the 
Julian date of the selected phase-zero point corresponding to 
the maximum light of the star {(p^-Q.Q), and P=375 days is the 
pulsation period of the star (Imai et al., 1997). 

Figure 3 compares the calibrated MIDI spectra for both pe- 
riods (2009/2011) with the ISO/SWS spectrum of RTVir ob- 
served 1996 July 20. The ISO/SWS absolute flux level of RTVir 
was scaled to the IRAS 12 and 25 fim photometry of the star. 
The MIDI flux level is -50% lower than that of ISO/SWS. Two 
possibilities could account for this difference. The first one is 
that the beam of an AT (~1.1" at 10 jim) is small compared to 
that of the IRAS satellite'*. Because of the relatively short dis- 
tance to RTVir (220 pc), and the large extent of its dust en- 
velope, some of the N-band flux might be lost outside the AT 
field-of-view. A second possibility that cannot be excluded is 
that the flux level has changed since the IRAS observation (per- 
formed in 1983), owing to cycle-to-cycle stellar variability. The 
AAVSO lightcurve shows an increase in the star visual ampli- 
tude by a factor of about 4 from the IRAS to the MIDI observa- 
tions. Moreover, there is a ~20% decrease in the MIDI flux from 
P83 (cycle 1) to P86 (cycle 3) (see Fig.3). Since the P83 and P86 
periods are close in phase, this loss in brightness can be due to a 
cycle-to-cycle variability of the object. 

3. Morphology of RTVir in tlie mid-infrared 

To constrain the morphology of the object, we used the software 
GEM-FIND (GEometrical Model Fitting for INterferometric 
Data; Klotz et al. 2012b). For the best reliability, the geometri- 
cal fitting requires the largest possible uv coverage (Klotz et al., 
2012b). To check for interferometric variability of the source 
(which may compromise the use of the full data set for the fit), 
we can compare the visibility measurements of data sets #3 and 
#6. The two data sets have very similar baseline lengths and po- 
sition angles, but were observed during two different cycles (see 
Table 1). The comparison between the two green upper curves 

* Keeping in mind that both 12 and 25yum IRAS photometry were 
used to scale the ISO spectrum. 



of Fig. 1 shows that no significant interferometric variability is 
found'', unlike what is found for the MIDI N-band flux measure- 
ments (see Fig. 3), which are dominated by the emission from 
the extended circumstellar regions. The fact that no significant 
cycle-to-cycle variability is found for the 30 m MIDI visibilities, 
which probe the inner atmospheric regions, reveals different mo- 
tions of the circumstellar environment from one spatial scale to 
the other. 

To account for the central star, we used a spherical UD 
model. We used a spherical Gaussian distribution to reproduce 
the optically thin dusty envelope of the star To account for 
the deviation from zero-differential phase in the 8-9 /urn wave- 
length region of the 128 m baseline measurement (see Fig. 2), 
we added an offset Dirac delta distribution representing an un- 
resolved component close to the star The diameter of the UD 
representing the central star and the FWHM of the Gaussian dis- 
tribution representing the dusty envelope are both wavelength 
dependent, in order to account for the opacity variation in the 
close stellar atmosphere and of the dust over the N-band range. 
Both the flux ratio of the envelope to the central star and the flux 
ratio of the unresolved component to the central star are also 
wavelength dependent, in order to account for the flux variation 
across the N-band for all three components. 

The fitting strategy is the same as described in Klotz et al. 
(2012b). Figure 4 presents the best-fitting geometrical model 
together with the calibrated MIDI visibility and differential 
phase measurements of RT Vir. Figure 5 shows the normalized 
intensity distribution of the best model at 10.5 //m. Figure 6 
gives the wavelength-dependent parameters of the best-fitting 
model. 

All visibility and 30 to 89 m basehne differential phase mea- 
surements are very well reproduced by the model. The similar 
trend of the model and the 128 m baseline differential phases 
in the 8-9 //m range provides one simple solution to the nature 
of the asymmetry*. It is assumed here that the asymmetry cor- 
responds to an unresolved component having a flux of about 
13 Jy from 8 to 8.5 fim (and about 2 Jy through the rest of the 
N-band), located at a projected distance of about nine stellar 
radii (12 AU away from the central star). However, the position 
of the unresolved component is not well constrained owing to 
the very limited uv-coverage. The result from the geometrical 
modeling shows that the unresolved component does not dis- 
rupt the sphericity of the close circumstellar environment (see 
Fig. 5), and could therefore be related to, e.g., an H2O or SiO 
maser clump (i.e. Richards et al. 201 1, 2012). 

The overall larger stellar diameter (compared to the 12.4 mas 
IC/K' -bands diameter) in the N-band (see in the upper left panel 
of Fig. 6) is probably caused by the presence of warm molecular 
layers of H2O and SiO (e.g. Ohnaka et al. 2005). Furthermore, 
the increase in the stellar diameter from 16 mas (1.3 stellar radii) 
to 21 mas (1.7 stellar radii) between 10.3 and 12.5 /im is pos- 
sibly correlated to the presence of AI2O3 (corundum), which is 
optically efficient around 11 //m (see Fig. 10), and has already 
been suggested to explain similar features seen in mid-infrared 
interferometric measurements of other comparable sources (e.g. 
Wittkowski et al. 2007; Zhao-Geisler et al. 2012). 



' With the exception of wavelengths longward of 12.3 /im that can be 
due to an underestimation of the uncertainties for data sets #3 and/or #6 
because of edge band noise within this wavelength range or because of 
the small baseline length difference (2 m) between both data sets. 

Keeping in mind that due to the lack of spatial coverage, more com- 
plex morphologies cannot be excluded. 
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Fig. 4. Best fit of the UD+Gauss+Dirac geometrical model (red solid 
lines) on the calibrated visibilities and differential phases of RTVir 
(black solid lines). 



The FWHM of the Gaussian distribution representing the 
dusty envelope increases from 16mas (1.3 stellar radii) to 
54 mas (4.3 stellar radii) between 8 and ~9.4yL(m. This increase 
is caused by a silicate dust opacity effect. The modeling also 
shows a large increase in the FWHM between 9.4 and 10.4/vm 
with a peak reaching 227 mas (18 stellar radii) around 9.8 yum. 
Such an increase is very likely related to an extended dusty 
circumstellar envelope composed of silicate grains. An extra 
shell of silicate dust also needs to be included in the dynamic 
modeling to explain the 30 and 32 m baseline visibility mea- 
surements (see Sect. 6). Finally, the FWHM shows a constant 
size of about 60 mas (5 stellar radii) between 10.4 and 12.5 yum. 

This first morphological interpretation already gives us some 
clues to the presence of aluminum and silicate dust in the close 
environment ( 1 .3 to 4.3 stellar radii) of RT Vir. Detailed theoret- 
ical modeling (Woitke, 2006) suggests that silicate grains form- 
ing on such spatial scales have to be iron-free not to evaporate 
from heating by the radiation field. 

The MIDI differential phases show the presence of an asym- 
metry in the stellar vicinity. However, since the asymmetry does 
not affect the overall sphericity of the close circumstellar en- 
vironment, this allows the use of the spherically symmetric dy- 
namical model to investigate the presence of aluminum and iron- 
free silicate dust. 



4. The dynamic model atmosphere of RTVir 

The dynamic model atmospheres (called DMAs hereafter) used 
are described in detail in Hofner (2008). They contain a con- 
sistent time-dependent description of grain growth for pure 
forsterite grains, and are based on radiative transfer in 3 1 9 wave- 
length points, taking size-dependent dust grain opacity into ac- 
count. Each dynamic model starts from a hydrostatic initial 
structure, and the effects of stellar pulsation are simulated by a 
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Fig. 5. Normalized intensity distribution of the best-fitting geometri- 
cal model at 10.5 fim. The offset white circle represents the unresolved 
component. North is up and east is left. 
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Fig. 6. Wavelength-dependent parameter values of the best-fitting ge- 
ometrical model. The upper left panel shows the stellar diameter. The 
K-band diameter of the star is superimposed as a dashed line. The up- 
per right panel represents the FWHM of the dusty envelope. The lower 
left panel represents the flux ratio of the dusty envelope over the central 
star, while the lower right panel shows the flux ratio of the unresolved 
component over the central star. 



variable inner boundary just below the stellar photosphere (pis- 
ton model). The dynamic code provides the temporally varying 
temperature and density distribution of the object. Based on the 
radial structure at given instances in time, atomic and molecular 
opacities are computed using the COMA code^ (Aringer et al., 
2009). COMA solves the ionization and chemical equilibrium 
equations at a given temperature and density combination for a 
set of atomic abundances assuming LTE. The calculated opac- 
ities are then introduced into a spherically symmetric radiative 



' This version of the COMA code does not take dust scattering into 
account. 



5 



Sacuto et al.: The wind of RT Virginis 



transfer code to calculate the emergent intensity distribution at 
any desired spectral resolution. This intensity distribution is used 
to derive the synthetic spectrum and visibility profile, which are 
compared to the observational data of the star 

4.1. Strategy 

The strategy of selecting the most appropriate DMA for RT Vir 
is similar to the one used by Sacuto et al. (201 1) for R Scl. 

1. Constrain the best set of stellar parameters, i.e. surface 
gravity, mass, effective temperature, and metallicity, which 
are used for the initial hydrostatic structure of the DMA. 
This is done by fitting hydrostatic MARCS models to 
the broadband photometric and ISO/SWS spectrometric 
measurements of the star 

2. Constrain the additional input parameters of the DMA, i.e., 
the amplitude of the piston and abundance of seed particles. 
This is done in two steps: 

a. Compute a sample of DMAs using the best-fit stellar 
parameters derived from the hydrostatic modeling. Different 
velocity amplitudes of the piston and abundances of seed 
particles are used to produce models with different mass-loss 
rates and wind velocities. Models producing mass-loss rates 
and wind velocities in agreement with the values found for 
RT Vir in the literature are selected for further analysis. 

b. Among the set of models selected in the previous step, 
constrain the best velocity amplitude of the piston and abun- 
dance of seed particles, by comparing the time-dependent 
synthetic dust-free spectra to the broadband photometric and 
ISO/SWS spectrometric measurements. 

These points are detailed in Sects. 4.2, 4.3.1, and 4.3.2, re- 
spectively. 

4.2. Determination of the steiiar parameters of RTVir 

The stellar parameters are determined by comparing hydrostatic 
stellar atmosphere models to the spectrophotometric data. The 
hydrostatic model that provided the best fit to the broadband 
B, V (simbad database) and J, H, and K photometric data from 
Kerschbaum & Hron (1994) and the ISO/SWS spectra of the 
star up to 5 /^m was chosen. This limit of 5 fim is fixed because 
the contributions from atmospheric molecular and dust emission 
become non-negligible at longer wavelengths for optically thin, 
warm, dusty objects like RT Vir (Kraemer et al., 2002), and that 
cannot be accounted for by the hydrostatic modeling. 

A least-squares minimization fitting was done using a grid 
of MARCS models (Gustafsson et al., 2008), which has a ba- 
sic sun-like chemical composition as listed by Grevesse et al. 
(2007). The eflFective temperature of the star was varied from 
2800 to 3200 K (ATeff=100 K) for metalHcity values of -0.25, 
0.0, and 0.25. The surface gravity log ^=-0.5 of the model was 
fixed according to a stellar mass M=l M©, the distance d=220pc 
deduced by Imai et al. (2003), and the angular diameter of the 
central star of 12.4 mas derived from K-band long-baseline in- 
terferometry by Monnier et al. (2004). Table 2 summarizes the 
stellar atmosphere characteristics of RTVir that was adopted, 
as well as deduced, from the best-fitting hydrostatic model. 
Figure 7 presents the best-fitting MARCS model, together with 
the spectrophotometric data. A discrepancy between the model 
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Fig. 7. Best-fitting hydrostatic MARCS atmospheric spectrum (green 
sohd line) and best-fitting dust-free spectrum of the DMA 1 (red solid 
line) at bolometric phase 0(,oi=O.59, together with the broadband photo- 
metric data (black squares) and ISO/SWS spectrum (black solid line) of 
RT Vir. The insets present the flux deviation of the hydrostatic MARCS 
spectrum (green solid line) and the flux deviation of the dust-free spec- 
trum of the DMA 1 (red solid line) from the ISO/SWS spectrum around 
3 fim and 4.25 //m, respectively. 



and the measurements is seen around 3 yum corresponding to the 
H2O and OH absorption bands, as well as around 4.25 fim cor- 
responding to the CO2 absorption band. The photospheric part 
of these bands provided by the hydrostatic model may not be 
enough to explain the observed features that could originate in 
the extended atmospheric environment. Finally, the temperature 
of Teff =2900 K leads to a luminosity of L=5500Lq for the star. 



Table 2. Stellar atmospheric parameters of RT Vir fixed and deduced 
from fitting a hydrostatic model to the spectrophotometric measure- 
ments. 



Method 


Parameter 


Value 




Stellar diameter: 9 (mas) 


12.4 


Fixed 


Stellar mass: M (Mq) 


1 




Surface gravity: log g 


-0.5 


Deduced from the fit 


Effective temperature: T^if (K) 
Metallicity: Z (Zq) 


2900 
1 



Now that the stellar parameters of RT Vir have been deter- 
mined, they are used for the initial hydrostatic structure to com- 
pute the corresponding DMAs having mass-loss rates and wind 
velocities in agreement with the ones found in the literature. 

4.3. Constraining fhe additional parameters for fhe DMA 

4.3.1. Dynamic model atmosphere sample 

In addition to the stellar parameters determined from the best 
fit of the hydrostatic model to the spectrophotometric measure- 
ments (see Sect. 4.2), the dynamical computation requires three 
other free parameters to be determined: (i) the abundance of 
seed particles per H atoms iigjnu, (ii) the velocity amplitude of 
the piston Aup, and (iii) the period of the piston Pmod- 

While the period of the piston is fixed to Pn,od=375 days, 
i.e., the pulsation period of the star (Imai et al., 1997), different 
values of ng,/nH and Aup are used, resulting in the development 
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of a wind, characterized by a mass-loss rate M and a wind 
velocity u. The range of values for Aup and ng^/nn are con- 
strained by the values of M and u derived from other studies: 1 
to 5x10"^ Moyr"' (Bowers & Johnston, 1994; Olofsson et al., 
2002; Imai et al., 2003) and 6.2 to llkms"' (Olofsson et al., 
2002; Gonzalez Delgado et al., 2003; Richards et al., 201 1). 

The parameters and resulting wind properties of the detailed 
DMAs are listed in Table 3. 



Table 3. DMAs parameters and resulting wind properties. 



DMA* 


Hgr/nn 


Aup 


<M> 


<u> 






[km s-'] 


[Moyr-i] 


[kms-'] 


1 


3x10^'" 


3 


1.9x10^' 


6.6 


2 


3x10-'^ 


3.5 


6.7x10^^ 


114 


3 


1x10-" 


4 


7.7x10^' 


7.1 


4 


1x10"'^ 


5 


1.4x10^* 


11.5 



The fundamental stellar parameters of the DMAs presented here are 
the ones listed in Table 2. 

'rigr/riH and AUp are inputs values of the model, while <M> and <u> 
are output values. 




-5.0 -4.5 -4.0 
log(grain radius) [cm] 



3.0 



Fig. 8. Absorption Qats (red solid line) and scattering 2.„,„ (blue dashed 
line) efficiencies of Mg2Si04 as functions of grain radius at /l=9.8yum. 
Data for the refractive index are taken from Jager et al. (2003). The 
efficiencies are calculated using the Mie theory for spherical particles 
(program BHMIE from Bohren & Huffman (1983) and modified by 
Draine, www.astro.princeton.edu/draine/scattering.html). 



4.3.2. Constraining tine DMA of RTVir 

To select the most appropriate DMA for RT Vir, among the ones 
given in Table 3, we compare the time-dependent spectra com- 
puted without dust opacities of each DMA 1/2/3/4 to the spec- 
trophotometric data below 8 /vm. This limit is fixed because of 
the ability of the dynamical models, in contrast to the hydro- 
static model, to self-consistently develop extended atmospheric 
molecular structure of H2O, SO2, and SiO typically emitting 
between 5 and 8 yum. Beyond 8 yum, the dust emission is non- 
negligible, so addition of dust opacities in the radiative transfer 
calculation is required (see Sect. 5). 

We investigated 46 snapshots from each model presented in 
Table 3, equidistant in phase^ and distributed over three consecu- 
tive pulsation cycles. The additional parameters, ngr/nH and Aup, 
were derived by fitting the average dust- free synthetic spectrum 
for each DMA (averaged over all phases and cycles) to the spec- 
trophotometric data below 8 /vm. Another test was done by find- 
ing the individual dust-free spectrum that presents the best agree- 
ment with the spectrophotometric data for a given bolometric 
phase. In both cases, DMA 1 gives the best agreement with the 
spectrophotometric data (see Fig. 7). Spectra of DMA 2/3/4 hav- 
ing a higher <M>, overestimate the flux compared to the spec- 
trophotometric measurements. While the best-fitting DMA 1 can 
reproduce the 3 yum band, in contrast to the hydrostatic model, it 
is not able to reproduce the measurements around 4.25 jum (CO2 
absorption band), as is also the case for the hydrostatic modeling 
(see Sect. 4.2). Because the contribution of extended molecular 
layers is accounted for by the DMA, the reason for the 4.25 yum 
discrepancy can be a deviation from chemical equilibrium in the 
low density outermost atmosphere (400-1 000 K), where the CO2 
feature is formed. In contrast to the hydrostatic MARCS models. 
Fig. 7 shows that the best-fitting DMA 1 is able to reproduce the 
6.5-8//m excess of H2O, SO2, and SiO emission from the ex- 
tended atmosphere. 



* The phase of the DMA corresponds to the bolometric phase (/iboi 
derived from its bolometric lightcurve (Nowotny et al., 2005b). 



5. Dust chemistry in tlie inner wind region of RTVir 

A major point of this study is to test the theoretical predictions 
of Hofner (2008), i.e., that winds of M-type AGB stars can be 
driven by photon scattering on iron-free silicate dust formed 
around two to three stellar radii. To do so, we compare the N- 
band spectro-interferometric measurements with the forsterite 
dust-rich spectrum and visibility profiles of the selected DMA 1 
(see Sect. 4.3.2). 

The predicted mean grain size of DMA 1 is <agi>=0.3 jum. 
Since the scattering efficiency of forsterite grains with sizes be- 
low 1 yum is lower than the absorption efficiency in the mid- 
infrared (see Fig. 8), this allows the use of COMA (with no scat- 
tering included) for the comparison between the DMA and the 
MIDI interferometric data. 

5.1. Iron-free silicate grains as wind drivers 

Using the temperature and density structures given by the dy- 
namical calculation, the intensity distribution can be computed 
within the N-band wavelength range at the MIDI-GRJSM spec- 
tral resolution (230) with the help of the radiative transfer code 
COMA for the 46 snapshots of DMA 1. Because the COMA 
code uses the same molecular and dust (forsterite) opacities, 
such as the ones included in the dynamical computation, the 
radiative transfer is done consistently. By Fourier transform of 
the intensity profiles for the different projected baselines of the 
observations (see Table 1), we can derive the corresponding dis- 
persed visibility profiles and compare them to the MIDI mea- 
surements of the star, following the same approach as described 
in Sacuto et al. (2011). We then select the best-fitting bolomet- 
ric phases of DMA 1 by least-squares minimization between the 
results of the model and the visibility measurements for periods 
2009 (cycle 1) and 201 1 (cycle 3). 

Figure 9 shows the best-fitting synthetic visibility profiles 
with forsterite of DMA 1, together with the MIDI mea- 
surements. We show, for comparison, the visibility profiles 
computed without dust opacities for the same bolometric phases 
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Fig. 9. Best-fitting synthetic visibility profiles with forsterite of DMA 1 
at phase 0boi=O.6O and 0(,oi=O.86 (red solid lines), together with the 
MIDI measurements of the star (black solid lines) for periods 2009 (cy- 
cle 1) and 201 1 (cycle 3), respectively. Red bars, next to the legend, rep- 
resent the average phase dispersion of the model. Red dotted lines repre- 
sent the corresponding dust-free visibility profiles at the same bolomet- 
ric phases. Unshaded areas represent the silicate dust's main effective 
domain. 



to demonstrate the impact of the silicate grains. 

The unshaded areas of Fig. 9 correspond to the wavelength 
region where silicate species show strong spectral features (see 
Fig. 10). Figure9 reveals that the fully self-consistent DMA 1, 
producing forsterite dust particles of sizes ~0.3 jim, is in better 
agreement with the interferometric measurements than the dust- 
free model, especially for the 59 and 63 m baselines, taking the 
average phase dispersion of the model into account shown by the 
bar next to the legend of each plot. 

This result provides evidence for the presence of silicate dust 
in the region around two to three stellar radii probed by the 
60 m baseline (considering an angular diameter of the star of 
12.4 mas). This gives observational support, in addition to what 
is provided by Norris et al. (2012), in favor of the theoretical 
predictions of Hofner (2008), i.e., that winds can be driven by 
photon scattering on iron-free silicate grains formed in the close 
vicinity of M-type AGB stars. A further test for the presence 
of silicate dust around two to three stellar radii is discussed in 
Sect. 5.3. 



5.2. Exploring extra dust opacities 

In the dynamical computation, the formation of only one 
type of dust, i.e. pure forsterite (Mg2Si04), is considered. 
However, it is likely that other types of dust can contribute 
to the enhancement of the opacity in the mid-infrared region. 
Figure 9 reveals that the synthetic visibilities are above the 
measurements longward of lOyum, indicating that material 
strongly emitting around 1 1 jum in the close stellar vicinity is 




10 11 
A, [urn] 



12 



13 



Fig. 10. Optical efficiency per grain radius Qots/ag,- as a function of 
wavelength in the mid-infrared region for a selection of dust species 
likely to exist in the circumstellar environments of M-type AGB stars. 
Refractive index data were taken from Begemann et al. (1997) (AI2O3), 
Zeidler et al. (2011) (TiOa, MgAl204), Palik (1985) (SiOa), Jager 
et al. (2003) (MgSiOj and Mg2Si64), and Dorschner et al. (1995) 
(MgFeSi04). 



missing in the radiative transfer computation. Figure 10 shows 
the efficiency per grain radius as a function of wavelength for 
a selection of oxygen-bearing dust species. From this plot, the 
only dust component that is optically efficient around 1 1 fim is 
AI2O3 (corundum) dust. Therefore, to improve the modeling 
of the visibility measurements longward of lOyum, we try to 
compensate for the lack of opacity of the model around 1 1 fim 
by adding alumina dust in the a-posteriori radiative transfer 
computation. 

Additional dust species can be included a-posteriori in the 
COMA radiative transfer code, as was done by Sacuto et al. 
(2011) for SiC dust in the case of RScl. The radiative trans- 
fer is done inconsistently in this case because for M-type AGB 
stars, the dynamical computation does not include other dust 
opacities than forsterite. However, because of the low abundance 
of Al, aluminum dust should not influence the process of wind 
driving in the dynamical computation (Bladh & Hofner, 2012a). 
Nevertheless, aluminum dust, like corundum, is likely to con- 
tribute to the opacity enhancement around 11 yum (see Fig. 10) 
leading to better agreement between the model and the mid- 
infrared visibility measurements. Finally, interferometric mea- 
surements in the mid-infrared have already pointed to the pres- 
ence of AI2O3 dust around two stellar radii in Mira and semi- 
regular M-type AGB stars (e.g. Wittkowski et al. 2007; Zhao- 
Geisler et al. 2012). We therefore decided to add corundum in 
the a-posteriori COMA radiative transfer computation, using the 
small-particle limit approximation. 

Refractive indices of corundum were adopted from 
Begemann et al. (1997), while the condensation temperature was 
fixed to rf^'^°' = 1400K (see e.g. Gail 2010). Since aluminum 
dust is not included in the dynamical computation, we assumed 
that corundum follows the same temperature-density structure 
as the gas. The only free parameter left is the degree of con- 
densation (/7A1.0,) of Al condensing into corundum. We varied 
77ai,03 from 10% to 100% (with Ari\uo^-lQ%) to minimize the 
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difference between the model and the visibiHties for all base- 
lines. Even if the addition of AI2O3 dust increases the opac- 
ity and leads to a decrease in the visibility amplitude longward 
of 9/vm, 100% of Al condensing into corundum is not suffi- 
cient to reach the level of the visibility measurements. However, 
because the overall shape of the synthetic visibility profile is 
similar to the measurements longward of 9 fim, we decided to 
increase the ratio above 100% assuming the creation of more 
dust than available Al atomic material, i.e., exceeding the so- 
lar abundance of Al. Varying 77ai,o, from 200% to 1000% (with 
A/7Ai2O, = 100%), we find the best agreement with the measure- 
ments for 77Ai2O3=500%. 

From the measurements, a clear decrease in the visibility 
profile is revealed above 12.5 ;um. Going back to Fig. 10, both 
MgAl204 (spinel) and Ti02 (titanium dioxide) show an increase 
in their optical efficiency longward of 12.5;um. Since TiOa has 
a lower optical efficiency than MgAl204 in the 12-13;um region 
(see Fig. 10), the contribution from the TiOa opacities will have 
a negligible impact in comparison to the MgAl204 opacities. 
Therefore, we decided to do the test by only adding spinel, in 
the same way as corundum, in the a-posteriori COMA radia- 
tive transfer computation, using the small-particle limit approx- 
imation. Refractive indices of spinel are adopted from Zeidler 
et al. (2011), while the condensation temperature is fixed to 
7'^^^'"°'' = 1 150K (see e.g. Gail 2010). Since aluminum dust 
is not included in the dynamical computation, we assume that 
spinel follows the same temperature-density structure as the gas. 
The degree of condensation (77MgAi204) is varied from 10% to 
100% (with A/7MgAi2O4 = 10%) to minimize the difference be- 
tween the model and the visibilities for all baselines. We found 
the best agreement for a full condensation (?7MgAiTO4 = 100%) of 
Al condensing into spinel. 

Figure 1 1 shows the comparison between the best-fitting 
synthetic visibility profiles with forsteri ten- aluminum of DMA 1 
and the MIDI visibility measurements. Visibility measurements 
for baselines equal to and larger than 59 m are reproduced well 
throughout the whole N-band range. 

Obviously, a degree of condensation above 100% (i.e. bind- 
ing up more of an element in dust than what is available) for 
Al or other elements is not a realistic solution, no matter how 
good the fit of the synthetic visibilities is to the observed ones. 
It should therefore be remembered that what is actually con- 
strained by the fit is the column density of dust material con- 
tributing opacity to the relevant wavelength region. The dust col- 
umn density is an integral containing the product of gas density 
times the degree of condensation. One possible solution to the 
problem is that the gas density in the circumstellar envelope in 
the close vicinity of RT Vir is actually higher than in our model. 
This would bring down the degrees of condensation by the cor- 
responding factor, while keeping the same dust column density, 
in order to fit the observations in the N band. Considering that 
the mass-loss rates taken from the literature (see Sect. 4.3.1) are 
uncertain by a factor of 3 or more, and are mean values span- 
ning more than several hundred years of mass loss, while RT Vir 
shows clear cycle-to-cycle variability and long-term trends (see 
Sect. 2.2), the current gas densities in the inner envelope are not 
well constrained by these values. 

While a higher gas density (corresponding to a higher mass- 
loss rate) may help solve problems in the dust-dominated spec- 
tral regions, it could, on the other hand, also affect the molecule- 
dominated regions of the SED, in particular in the near-infrared. 
However, as illustrated by the current problems reproducing the 
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Fig. 11. Best-fitting synthetic visibility profiles with 
forsterite -I- aluminum of DMA 1 at phase (^boi=0.60 and (^boi=0.86 
(red solid lines) together with the MIDI measurements of the star 
(black solid lines) for periods 2009 (cycle 1) and 2011 (cycle 3), 
respectively. Red bars, next to the legend, represent the average phase 
dispersion of the model. Red dotted lines represent the synthetic 
visibility profiles with aluminum dust only of DMA 1 at the same 
bolometric phases. 



CO2 feature (see Sect. 4.3.2), questions remain regarding the va- 
lidity of chemical equilibrium in the gas phase (e.g. Cherchneff 
2006) and other model assumptions which could affect the 
molecular spectra. Therefore a discussion of this problem is be- 
yond the scope of the current paper. 

5.3. Identification of silicate grains in the 2 to 3 stellar radii 
region 

To confirm the presence of silicate dust in the two to three stel- 
lar radii region, the following test was performed. The COMA 
radiative transfer code was run only including aluminum dust 
(with no silicate dust included) for the 46 snapshots of DMA 1 . 
Corundum and spinel dust grains were added in the same way as 
in the Sect. 5.2. Refractive indices of corundum and spinel were 
adopted from Begemann et al. (1997) and Zeidler et al. (2011), 
respectively, while the condensation temperature was fixed to 
T=1400K for corundum and T=l 150K for spinel. The only free 
parameter left is the degree of condensation of Al condensing 
into corundum and spinel. For that we varied the degree of con- 
densation from 100% to 1000% (with steps of 100%) to mini- 
mize the difference between the model and the visibilities for all 
baselines. 

Figure 1 1 presents the best-fitting synthetic visibility pro- 
files with aluminum dust only of the DMA 1 and the visibil- 
ity measurements of RT Vir. We found the best agreement with 
the measurements for the same degree of condensation of Al 
(77AI2O3 =500% and /7Mi;AiiO4 = 100%) as the model that includes 
both forterite and aluminum dust (see Sect. 5.2). 
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Figure 1 1 shows a clear discrepancy (taking the average 
phase dispersion of the model into account) between the syn- 
thetic visibility profiles only with aluminum dust and the mea- 
surements around 9.8 //m for the two 59 and 63 m baselines 
measurements''. This reveals a lack of opacity within the re- 
gion around two to three stellar radii probed by these base- 
lines at 9.8 fim. Among the data shown in Fig. 10, only the sili- 
cate species show strong spectral features around 9.8 fim. From 
Fig. 10, it is difficult to distinguish between iron-rich and iron- 
free silicate dust in this wavelength region. However, detailed 
theoretical modeling (Woitke, 2006) suggests that silicate grains 
forming within the region around two to three stellar radii have 
to be iron-free in order not to evaporate due to heating by the ra- 
diation field. Figure 1 1 shows that the model that includes both 
forsterite and aluminum dust (with the same degree of condensa- 
tion of Al) can indeed reproduce the 9.8 ;um visibility measure- 
ments for the two 59 and 63 m baselines, giving a clear confir- 
mation for the presence of silicate in the two to three stellar radii 
spatial region. 

The model including only aluminum dust (see Fig. 11) also 
reveals that silicate dust has a very slight impact on spatial scales 
around two stellar radii (probed by the 89 m baseline) and almost 
no influence on spatial scales around 1 .5 stellar radii (probed by 
the 128 m baseline), taking the average phase dispersion of the 
model into account. This is reasonable assuming that the heat- 
ing by the radiation field within this region certainly leads to 
the evaporation of silicate dust but not to the evaporation of alu- 
minum dust. 

This result not only confirms the presence of silicate dust in 
the region around two to three stellar radii, but also indicates that 
the column density of forsterite dust, provided by the DMA 1 , 
is able to reproduce the spatial extent of the inner wind region 
within the 8-10yL(m spectral range. 



6. Exploring geometrical issues 

Figure 1 1 reveals a lack of opacity of the model in compari- 
son to the smallest 30 m baseline visibility measurements. Such 
an opacity enhancement was also found around 9.8/vm by the 
N-band geometrical modeling of RT Vir (see Sect. 3 and upper- 
right plot of Fig. 6). One possibility for overcoming this discrep- 
ancy is to artificially add an extra shell of dust material. 

6.1. Addition of an extra she// of forsterite 

We decided to include an extra shell of forsterite dust, in addition 
to our previous forsterite plus aluminum model (see Sect. 5. 2), in 
order to reproduce the 30 m baseline visibility data. To do that, 
we varied the position of the inner edge of the shell between four 
and ten stellar radii, corresponding to the spatial region probed 
by the 30 m baseline. The abundance of forsterite is also var- 
ied under the assumption that 20, 40, 60, and 80% Si has con- 
densed into forsterite. The best fit is obtained for the location 
of the extra shell around four to six stellar radii away depend- 
ing on the phase. The best forsterite abundance ratio is found to 
be /7MgTSiO4=60%, while the abundance of aluminum grains are 
found to be 77ai,o,=400% and /7MgAi,O4=80%, keeping the same 
condensation temperatures of 1400 K and 1 150 K, respectively. 



' Increasing the degree of condensation of Al condensing into corun- 
dum further makes the visibiUty profile of the model too steep in com- 
parison with the measurements in the 10-13 fim region without correct- 
ing the discrepancy around 9.8 fim. 




wavelength (/im) wavelength {/im) 

Fig. 12. Best-fitting synthetic visibility profiles with 
forsterite -I- aluminum and extra shell of forsterite of DMA 1 at 
phase 0boi=O.6O and 0boi=O.86 (red solid lines), together with the MIDI 
measurements of the star (black solid lines) for periods 2009 (cycle 1) 
and 201 1 (cycle 3), respectively. Red bars, next to the legend, represent 
the average phase dispersion of the model. 



Results are shown in Figs. 12 and 13. Visibilities are now 
reproduced relatively well for all baselines, even if some dis- 
crepancies with the measurements still remain. The coiTespond- 
ing spectrum also shows a relatively good agreement with the 
measurements over the whole spectral range, especially reveal- 
ing that the extra opacities added in the radiative transfer com- 
putation do not lead to any discrepancies in the near-infrared 
wavelength range. 

Figure. 14 presents the corresponding intensity profiles at 
10 fim, with forsterite plus aluminum dust and the extra shell of 
forsterite, of DMA 1 for phases 0boi=O.6O and 0boi=O.86. The ex- 
tra shell is distinguishable from the horn-shaped profile around 
four stellar radii at phase (/>boi=0.60, and around six stellar radii 
at phase 0boi=O.86. It is interesting to note that the position of 
the shell reproducing the cycle 1 interferometric measurements 
is around four stellar radii, while the position of the shell repro- 
ducing the cycle 3 interferometric measurements is around six 
stellar radii. This can be interpreted as shell having moved for- 
ward between the two cycles with a velocity of about 7.5 km s"' 
(in good agreement with the time-average wind velocity of the 
best DMA 1, <u>=6.6 kms"'). 

6.2. Possible explanations of the results 

Two scenarios, both related to the geometry of the circumstellar 
environment, can explain the need to add an extra dust shell to 
reproduce the interferometric measurements. 

- Radial structure 

Because of the variability of the object, radial structure 
changes with time. The visual magnitudes of the star reveal 
an amplitude increase in its lightcurve over time. Moreover, 
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Fig. 13. Best-fitting synthetic spectrum with forsterite+aluminum and 
extra shell of forsterite of DMA 1 dispersed in phases (red solid line 
+ bars), together with the broadband photometric data (black squares), 
ISO spectrum (black solid line), and MIDI spectra of RT Vir at phase 
0v=O. 15 (Cycle 1) (blue eiTor bars) and phase (^v~0.0 (Cycle 3) (light 
green error bars). The best-fitting hydrostatic MARCS atmospheric 
spectrum (green solid line) is added too. 
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Fig. 14. Corresponding synthetic intensity profiles at lOyum, with 
forsterite-l-aluminum and the extra shell of forsterite, of DMA I at 
phases (f)f,a\=0-60 (black solid line) and 0boi=O.86 (red dashed line). 



the comparison between the MIDI and ISO/SWS spec- 
tra shows a cycle-to-cycle variability (see Sect. 2.2). These 
changes can be related to a long-term variability of the star. 
Since we are only using a limited time-sequence of the 
DMA, we could be missing a period for which the density 
structure of the DMA gives better agreement with the inter- 
ferometric measurements, without the need for extra silicate 
opacities. 

- Deviation from spherical symmetry 

Differential phase interferometric measurements reveal the 
presence of an asymmetry in the atmosphere of RT Vir (see 
Sect. 3). This is likely related to the presence of one or sev- 
eral clumpy structures in the close vicinity of the star, also 
suggested from water masers observations (Richards et al., 
2011). Three-dimensional radiation hydrodynamics simula- 
tions of the convective interior and stellar atmosphere sug- 
gest there are non radial structures in the dust shells of 



AGB stars (Freytag & Hofner, 2008) that might explain these 
clumpy structures and could potentially affect the dynamical 
process of wind formation. 

7. Conclusion and perspectives 

This work is a first attempt to understand the wind mechanism 
of M-type AGB objects by comparing photometric, spectro- 
metric, and interferometric measurements with state-of-the-art, 
self-consistent, dust-driven wind models. 

Comparison between spectro-interferometric data and a self- 
consistent dust-driven wind model reveals that silicate dust has 
to be present in the region between two to three stellar radii 
to reproduce the 59 and 63 m baseline visibility measurements 
around 9.8 jim. Detailed theoretical modeling (Woitke, 2006) 
suggests that silicate grains forming on such spatial scales have 
to be iron-free in order not to evaporate owing to heating by the 
radiation field. 

This result also indicates that the column density of forsterite 
dust, provided by the DMA, is able to reproduce the spatial ex- 
tent of the inner wind region in the 8-10;um spectral range. This 
gives observational evidence, in addition to what is provided by 
Norris et al. (2012), in favor of the theoretical predictions of 
Hofner (2008), i.e., that winds can be driven by photon scatter- 
ing on iron-free silicate grains formed in the close vicinity of an 
M-type AGB star 

However, the synthetic visibilities resulting from the model 
are above the measurements longward of 10 fim, indicating that 
material strongly emitting around 1 1 fim in the close stellar 
vicinity is missing in the radiative transfer computation. Ad-hoc 
addition of corundum and spinel dust in the a-posteriori radiative 
transfer can compensate for the lack of opacity, giving better 
agreement between the model and the 10-13//m interferometric 
measurements for baselines longer than 59 m. 

The necessity to increase the degree of condensation of Al 
material to a value about five times higher than the solar abun- 
dance, in order to obtain a good fit to the MIDI measurements, 
may indicate that the circumstellar gas density in the close vicin- 
ity of RT Vir is actually higher than in our model. The mass-loss 
rates used to constrain the dynamic model atmosphere of RT Vir 
are uncertain by a factor of 3, or more, and are mean values span- 
ning more than several hundred years (e.g. Ramstedt et al. 2008). 
Considering that RTVir shows clear cycle-to-cycle variability 
and long-term trends, the current gas densities in the inner en- 
velope are not well constrained by these average values. On the 
other hand, a higher gas density could also affect the molecule- 
dominated regions of the SED, in particular in the near-infrared. 
However, as illustrated by the inability of the model to repro- 
duce the CO2 feature, questions remain regarding the validity of 
chemical equilibrium in the gas phase and other model assump- 
tions that could affect the molecular spectra. 

The need for an extra silicate dust shell to reach the right 
level of the 30 m baseline visibility measurements is still 
puzzling, but can be related to the geometry of the circumstellar 
environment in various ways. Among them, the long-term 
variability of the star and the presence of clumpy structures in 
the stellar vicinity are the most likely ones. To overcome these 
problems, a larger statistical study and progress in advanced 
self-consistent 3D modeling are required. 
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